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A  reaction-diffusion  model  suggested  by  Yakhno  as  a  model  for  the  centers 
of  periodic  activity  often  seen  in  excitable  media  such  as  cardiac  tissue  and 
special  chemical  reagents  is  subjected  to  multiple-scale  analysis.  A  fairly 
complete  description  of  the  process  (which  involves  the  generation  of  wave 
fronts)  is  obtained,  but  it  is  shown  to  be  structurally  unstable  as  it  depends 
on  symmetries  in  the  basic  equations.  The  allowed  amount  of  deviation  from 
symmetry  is  estimated. 
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SIGNIFICANCE  AND  EXPLANATION 


Leading  centers  are  locations  often  observed  in  excitable  media  such  as 
cardiac  tissue  and  special  chemical  reagents,  which  apparently  serve  as  the 
source  of  wavelike  activity  such  as  rotating  spiral  configurations  in  the 
chemical  concentration  or  physico-chemical  state.  A  possible  model  to  explain 
this  is  subjected  to  asymptotic  analysis. 
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ON  YAXHNO'S  MODEL  FOR  A  LE  '  .NG  CENTER 


Paul  C.  Fife 

1.  Introduction 

In  [1] ,  Yakhno  presented  a  rather  simple  class  of  reaction-diffusion  systems  which 
exhibit  a  phenomenon  he  called  "division"  of  arrested  wave  fronts,  resulting  in  the 
successive  generation  of  diverging  pairs  of  fronts  from  a  central  point.  These  divergino 
pairs  form  an  antisymmetric  composite  profile,  so  differ  from  such  front  generation 
processes  studied,  for  example,  in  [4],  15) ,  and  other  places.  This  mechanism  was  proposed 
as  a  possible  model  for  the  center  of  periodic  activity  sometimes  seen  in  such  excitable 
media  as  heart  muscle  (fibrillation)  and  the  Belousov-Zhabotinskii  reagents  (especially 
spiral  waves),  and  sometimes  called  leading  centers.  Due  to  the  great  interest  in  such 
centers,  it  is  appropriate  to  examine  Yakhno' s  mechanism  closely. 

The  model  was  further  discussed  in  (2)  and  [3] ,  where  the  results  of  numerical 
simulation  were  given.  It  was  brought  out  in  those  papers  that  when  the  symmetry 
properties  present  in  Yakhno' s  initial  example  are  violated,  the  mechanism  he  described 
might  not  occur,  due  to  "running  away"  of  the  wave  fron  jfore  it  has  a  chance  to  divide. 

The  purpose  of  the  present  paper  is  to  give  a  qualitative  analytical  treatment 
(without  computer  simulation)  of  Yakhno' s  model,  more  detailed  than  that  in  the  above 
references.  In  particular,  scaling  and  order-of-magnitude  arguments,  which  were  not  fully 
exploited  there,  are  used  throughout  the  paper.  In  our  description  of  the  division 
process,  the  known  global  stability  properties  of  fronts  and  similar  structures  are  an 
important  part  of  the  argument.  Finally,  the  question  of  whether  an  arrested  front  (which 
serves  as  initial  condition  in  the  processes  analyzed  in  11-3))  can  arise  in  the  medium 
from  initial  data  of  some  degree  of  generality  is  addressed  in  sec.  4. 
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The  general  conclusion  is  that  Yakhno's  mechanism  is  rather  structurally  unstable. 

This  appears  to  apree  with  the  tenor  of  the  discussion  in  [2]  and  (3).  The  equations 
originally  proposed  and  analyzed  by  Yakhno  have  some  symmetry  properties;  they  also  have 
two  natural  time  scales  (a  fast  and  a  slow  one).  If  symmetry  is  now  violated  in  the 
equations,  its  violation  has  to  be  by  a  small  amount,  in  some  sense  comparable  to  e  (the 
ratio  of  time  units  in  the  two  scales),  in  order  for  the  first  "division"  to  take  place. 

The  cumulative  effect  of  these  deviations  would  make  subsequent  divisions  even  more 
unlikely.  But  even  beyond  this,  our  analysis  of  the  effect  of  different  Initial  data  is 
also  discouraging.  Because  of  the  relaxational  nature  of  the  basic  eouations,  initial  data 
will  typically  rapidly  equilibrate  to  the  shape  of  a  wave  front,  whose  trigger  velocity  is 
a  function  of  the  value  of  the  slow  variable  v  at  the  front.  If  v  eguals  a  special 
value  vg  at  the  front,  we  have  the  arrested  front  described  in  the  cited  references.  But 
for  v  y  Vg  there,  we  show  in  sec.  4  that  regardless  of  the  distribution  of  v  away  from 
the  front,  the  latter  cannot  later  evolve  into  an  arrested  front,  which  is  necessary  for 
division  to  occur. 

The  effect  of  the  size  of  the  spatial  domain  was  considered  in  [3],  but  ignored  here; 
for  simplicity  we  take  it  to  be  infinite. 

Calculations  reported  in  [21  suggest  that  when  the  perturbation  from  symmetry  in  f 
is  large  enough,  conditions  suitable  for  division  may  be  restored.  No  explanation  has  been 
advanced  for  such  a  phenomenon. 

Finally,  despite  the  structural  instability  of  Yakhno's  mechanism,  it  is  entirely 
possible  that  a  more  complex  model,  in  which  a  control  process  is  adjoined  to  the  basic 
equations,  may  produce  the  required  stability. 

I  am  indebted  to  4rt  Winfree  for  stimulating  my  interest  in  this  mechanism. 


The  model  in  the  presence  of  symmetr 


The  model  Yakhno  proposed  is  a  pair  of  reaction-diffusion  equations 
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eu  »  a  u  +  f(u,v) 

t  XX 


v  -  g(u,v)  , 


(1) 

(2) 


where  0  <  e  <<  1,  f  has  the  nullcline  shown  in  Figure  1,  and  g  is  discontinuous  in 
u  at  u0  (this  was  not  assumed  in  later  analysis),  changing  sign  there  from  (-)  to 
(+).  For  the  qualitative  results  in  this  section,  there  is  no  loss  of  generality  in  taking 


g(u) 


for  u  >  u0 
for  u  <  u0 


(3) 


u 


Wo  suppose  the  spatial  domain  to  be  the  entire  x-axis,  -«•  <  x  <  “,  so  the  spatial 
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variable  may  be  scaled  at  will:  this  only  affects  the  value  of  the  coefficient  a  and  th 
form  of  the  initial  data  (if  it  is  an  initial  value  problem).  The  scaling  which  most 
readily  provides  a  clear  description  of  the  process  is  that  for  which  a  ■  e,  and  we  s*-al 
accordingly  take  it  to  be  so. 

Besides  the  symmetry  in  (3),  we  also  assume  in  this  section  that  the  function  f  is 
antisymmetric  with  respect  to  the  point  (uq.Vq): 

f(Ug  -  4u,  vQ  -  6v)  »  “f(Ug  +  6u,  Vg  +  6v)  .  '4! 

Finally,  we  suppose  that  the  initial  data  also  have  symmetry  properties: 

u(x,0 )  -  •flx/e.)  <  \  ) 

v(x,0)  «  vQ  ,  if  ' 

where  r(£)  is  bounded  with  bounded  derivatives, 

\  U'<5>  +  v'  (-5)1  ■=  uQ 


rind  £(*(£)  "  UQ  )  >  0  for  ^  ¥  0. 

We  shall  describe  the  various  stages  in  the  evolution  of  the  solution  u(x,t),  v  ( 
Stage  1 :  Equilibration  of  the  initial  data  to  a  stationary  wave  front  (duration 
0(c)),  Changing  to  variables  r  *=  t/e ,  %  *  x/e  ,  we  may  write  (1)  (with  a  =  e)  t<> 
lowest  order  as 


u  *  u»r  ♦  f(u,v_) 

T  \ r.  0 
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Here  we  have  used  the  fact  that  vfc  is  bounded,  so  in  .he  T  time  scale,  v  does  not 
deviate  appreciat)'/  from  its  original  value  vq.  Since  (by  symmetry) 


fw 

^h_(v0)  ftu'Vdu  “  0  ' 


we  know  that  (8)  has  a  globally  stable  stationary  wave  front  solution  u  =  U(C)  satisfying 
U(4“)  “  ht*vo*'  “  ug"  Furthermore,  by  [6],  the  solution  of  (8)  with 

u(£,0)  «  v>(5)  approaches  U  uniformly  in  £,  as  x  ♦  ».  So  in  a  T-time  interval  of 
0(1)  (t-interval  of  0(e)),  u(£,x)  will  be  attracted  to  a  neighborhood  of  U(£). 

Stage  2 :  The  evolution  of  v  and  the  concomitant  changes  in  the  profile  of 
u( duration  0(1)).  The  wave  profile  to  which  u  has  evolved  is  such  that  u  >  Uq  exactly 
when  x  >  0,  and  this  remains  so,  because  of  symmetry.  Therefore  according  to  (2),  (3), 
we  have 


v(  x ,  t ) 


x  >  0 
x  <  0  . 


(9) 


For  each  fixed  time  t  during  this  stage,  v  is  a  step  function  in  x. 

Again  because  of  the  symmetries,  it  will  be  clear  that  v  -  vQ  and  u  -  ug  will 
continue  to  be  odd  functions  of  x  for  all  later  time.  For  this  reason,  it  suffices  to 
consider  the  problem  only  on  the  half-line  -•»  <  x  <  0,  with  the  boundary  condition 
u(0,t)  =  Uq.  Of  course  no  such  boundary  condition  can  be  applied  to  v,  because  v  is 
discontinuous  at  x  =  0;  and  no  boundary  condition  is  needed  in  any  event,  because  the 
differential  equation  for  v  involves  no  x-derivatives.  The  following  result  is  relevant 
to  the  study  of  this  problem.  The  stability  spoken  of  in  the  lemma  is  C°  -  stability  o- 
the  half-line  (-“.O). 

Lemma.  Let  f(u)  be  defined  and  continuously  differentiable  on  the  interval 
(0,1),  with  f ( 0 )  =  0  and  f’(0)  <  0.  Let  8  e  (0,1).  The  problem 
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We  apply  the  lemma  to  the  problem  of  finding  stable  stationary  solutions  u  =  U(£;v), 


£  =  i  of  (1)  with  a  =  e,  v  =  const,  on  the  half-line  (-“,0),  satisfying 
U(-";v)  =  h_(v),  U(0;v)  =  ug* 

Referring  to  Figure  1,  we  see  that  there  exists  such  a  solution  if  and  only  if 
v  >  v_-j,  the  latter  defined  by  the  condition 


0  . 


In  stage  1,  as  we  have  seen,  the  u-profile  is  rapidly  attracted  to  U(£;v0>.  During 
the  present  stage,  v  is  steadily  evolving  according  to  (9).  At  the  same  time,  u(x,t) 
continuously  adjusts  to  the  stable  profile  U(£;  v  -  t),  the  characteristic  time  variable 
associated  with  the  adjusting  process  being  the  fast  time  t.  In  short,  during  stage  2, 
for  x  <  0, 


v  (  X ,  t )  ~  v0  -  t  , 
u(x,t)  ~  IH-^r  vQ  -  t) 


(10) 


This  stage  continues  only  so  long  as  vQ  -  t  >  v_i-  For  time  t  >  v0  -v_v  according 
to  the  lemma,  the  equation  (1)  no  longer  has  a  stationary  solution  satisfying 
IK-”)  =  h  ( v)  (the  two  lower  curves  in  Figure  2  annihilate  themselves). 

Stage  3:  The  division  of  the  front  (duration  a  time  interval  <<  1).  For  v  <  v_ 1 , 

the  only  stationary  attractor  for  (1)  satisfying  the  boundary  condition  at  x  =  0  is  the 

monotone  decreasing  solution  u  =  U(£;v)  satisfying  (J(£;v)  =  h+(v)  (the  upper  curve  in 

Figure  2),  and  the  u-profile  is  now  attracted  to  that  solution.  This  abrupt  flip  from  one 
stable  state  to  another  is  reminiscent  of  a  variety  of  such  events  known  to  occur  in 
chemical  engineering  and  chemical  physics,  when  multiple  stable  steady  states  occur 
( ' F ,  q,  11],  for  example).  When  diffusion  is  present,  the  "flip”  characteristically 
generates  wave  fronts.  In  the  present  case,  consider  a  profile  for  which  n(-«°,v)  =  h  (v) 
which  is  in  the  domain  of  attraction  of  U(£,v).  For  large  negative  £ ,  u  takes  on 
va’-ins  near  h_(v),  which  is  a  stable  zero  for  the  function  f(u,v).  The  function  u 


f  or 


cannot  rapidly  depart  from  that  neighborhood  of  h_(v).  However,  there  is  a  mechanism 
the  eventual  accomplishment  of  that  goal.  If  v  were  held  constant  (rather  than  bein.o 
discontinuous  at  x  =  0),  there  would  exist  a  wave  front  solution  u(-,T)  =  W  ( '  -  c(v)T/, 
c  <  0,  satisfying 


ut  =  u^  +  f ( u , v ) ,  u(±°»,T)  =  h±(v) 

(the  negativity  of  c  follows  from  the  fact  that  the  state  h+  is  dominant,  i.e. 
h+  ( v ) 

f ( u,v) du  >  0  in  the  terminology  of  [9]  (see  also  [10])).  This  wave  front  is  very 

'  h  ( v) 

stable  (6).  What  we  will  see  in  Stage  3  is  the  initiation  of  a  structure  which  is  not 
precisely  this  wave  front,  but  (in  the  x  —  t  coordinate  system  again)  is  a  sharp  front, 
propagating  to  the  left,  attached  by  a  boundary  layer  to  the  fixed  value  u„  at  x  =  0 , 
as  shown  in  Figure  3.  It  will  also  be  seen  below  that  the  velocity  of  propagation  will  no 
be  constant,  as  v  is  not  constant. 


h  (v) 


> 


Figure  3 


It  is  intuitively  clear  (and  no  doubt  provable;  that  ...is  is  a  stable  configuration,  and 
is  the  mode  to  .  ich  the  solution  "switches"  when  v  become  less  than  v_.|.  Convergence 
of  solutions  to  similar  structures,  with  v  constant,  was  studied  in  (13; .  Stage  3  is  t 
initial  stage,  in  which  the  left-facing  front  is  first  formed  (the  dotted  line  in  "igure 
3). 

Stage  4;  The  propagation  of  the  divided  front  (time  scale  0(1)).  The  velocity  of 
the  propagating  wave  shown  in  Figure  3  is  seen,  by  dimensional  analysis,  to  be  0(1)  in 
the  (x  -  t)  system,  as  well  as  in  the  (£  -  t)  system.  In  fact,  it  is  clearly  0(1) 
the  latter  system,  and  so  the  traveling  wave  would  behave  like  a  function  of 

5-cx=— -c— =  —  —  (c  =  0(1)).  The  width  of  the  front  is  0(e)  since  the 
e  e  e 

characteristic  spatial  variable  is  C  = 

We  have  now  (Figure  3)  the  picture  of  a  sharp  front  moving  to  the  left  with  velocity 
0(11,  and  at  the  same  time  v  is  evolving  at  a  rate  0(1).  In  fact  by  (3),  v  is 
increasing  at  a  constant  rate  on  the  interval  between  the  front  and  the  origin,  and 
decreasing  on  the  left  of  the  front.  Letting  y(t)  be  the  position  of  the  front  (where 
u  =  Uq ,  say)  and  t1  =  vQ  -  v_,  the  time  at  which  sta-  -  3  occurs,  we  can  calculate  the 
approximat  position  of  the  front,  and  the  distribution  of  v,  as  functions  of  t  as 
follows.  To  the  left  of  the  front,  v  is  steadily  decreasing,  so 

v(x,t)  ~  v  -  (t-t^,  x  <  y(t)  .  (11) 

The  front's  velocity  cF  is  determined  by  this  value  of  v,  since  the  velocity  of  the 
front  solution  of  (8)  can  depend  only  on  v  in  that  equation  (10],  which  is  the 
(relatively  well  defined)  value  of  v  at  the  front.  As  in  (4) ,  we  therefore  have 

cF(  t)  =  c( v_1  -  ( t-t 1 ) ) 

Note  that  c  <  0  and  the  front  accelerates,  since  |c|  is  a  decreasing  function 
(see  Appendix  in  [5)).  Integrating,  we  have 


y  =  y  ( t ) 


c._  ( s)  ds 


Inverting  this  function  yields  T(xl#  the  time  of  arrival  of  t ho  front  at  position  x. 
Finally#  the  function  v  is  given  by 


v(x,  t)  =  v  -  (T(  x )  -  )  +  ( t-T  (  x ) )  #  y{  t )  <  x  <  0  .  M2) 

Stage  5:  The  vanishing  of  the  moving  front  (duration  <<  1).  Eventually#  the  value 
of  v  at  the  front  will  reach  v,  the  minimal  value  on  the  nullcurve  in  Figure  1.  For 
v  <  v^  h_(v)  no  longer  exists#  and  the  front  can  no  longer  exist.  When  this  happens#  u 
will  approach  the  value  h+(v),  uniformly  for  all  x  <  $  <  0#  effectively  annihilating 
the  propagating  front.  The  characteristic  time  scale  for  this  process  is  <<  1.  At  this 

point#  the  distribution  of  v  is  given  by  (12)  with  t  «  t2,  the  time  of  Stage  b.  Of 

course  v(x#t2)  =  v  for  x  <  y(t2).  Notice  that  v(x#t2)  attains  a  maximum  of 
2 v_ ■)  -  _v  at  x  -  0. 

Stage  6.  The  reversed  stationary  front  (duration  0(1)).  Now  u  >  uQ  everywhere 
for  x  <  0#  and  so  v  increases  with  uniform  velocity  1.  Just  as  in  Stage  2,  the  u- 
profile  continuously  adjusts  to  the  increasino  value  of  v.  In  the  boundary  layer  just  to 
the  left  of  the  origin#  the  sharp  variation  of  u  is#  in  the  stretched  coordinates,  a 

function  decreasing  from  ^'v)  to  up.  Hence  the  value  of  v  there  is 

v'x-0#t)  -  2v_,j-  v  +  bsitsid'*  the  boundary  layer,  of  course#  f(u#v)  -  n , 

so  u(x,t)  ~  h(vfx#t))#  ar.  i  v  is  determined  as  indicated  above. 

"""is  stage  continues  ..util  vft)  reaches  the  value  v1 . 

Stage  7:  Second  front  division  'duration  <<  1).  This  division  produces  a  downjuno 
front  moving  to  the  loft,  analogous  to  the  up  jump  produced  in  Stage  3. 

Remark :  We  have  assumed  that  Stage  s  occurs  before  Staae  ,  but  in  some  case.;  the 
order  will  be  reversed.  In  such  cases#  there  will  simultaneously  exist  an  up jump  front  and 
a  in vnjimn  moving  to  the  left,  bof^r*-  the  front  is  eliminated.  If  tne  function  g  in  (2) 
is  altered  so  that  ;t  Jer>er,ds  o»>  /  as  wolJ,  and  is  such  that  n  =  0  for  some  »ra]n« 


-in- 


v*  e  (v,v  )  and  u  <  uQ,  then  Stage  5  will  never  he  reached,  and  an  infinite  train  of 
fronts  will  accumulate. 

Of  course  to  the  right  of  the  origin,  a  succession  of  fronts  (of  the  opposite  shape) 
are  generated  in  the  antisymmetric  fashion.  Thus  at  all  times,  the  picture  is  such  that  a 
stationary  front,  facing  either  left  or  right,  exists  at  the  origin,  and  in  addition  fronts 
may  exist  which  diverge  from  the  origin  in  opposite  directions.  The  whole  process  is  shown 
schematically  in  Figure  4.  The  profiles  shown  there  are  the  u-profiles  at  the  end  of  the 
various  successive  stages.  The  numbers  indicate  the  stages. 
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We  now  consider  whether  the  phenomenon  of  front  division,  as  described  in  sec.  2, 


'-ar. 


i.ul'i  occur  when  the  various  special  assumptions  made  there  are  modified.  For  tnis 
pcijose,  we  first  ask  what  features  of  the  process  are  the  essential  ones  for  front 
ii'isicn,  and  then  whether  these  features  are  preserved  or  lost  when  the  symmetry 
conditions  are  relaxed. 

The  most  essential  feature  seems  to  be  that  the  function  v(x,t)  be  capable  of 
achieving  a  spatial  variation  of  the  order  of  magnitude  1  within  the  confines  of  a  single 
wave  front.  In  the  symmetric  case,  v  has  a  jump  discontinuity  of  this  magnitude  exactly 
at  the  center  of  the  front,  which  has  width  0(e),  and  so  that  property  holds.  Of  course 
it  is  not  necessary  for  v  to  be  discontinuous,  but  it  should  at  least  vary  by  such  an 
amount  in  an  interval  of  width  0(e)  centered  at  the  front.  In  fact,  if  it  does  not,  the 
to  first  approximation,  v  may  be  considered  constant  through  the  front,  and  the  front 
would  represent  a  horizontal  transition  in  Figure  1  from  the  branch  h_  to  h+.  Such 
•ronts  always  exist,  whereas  it  is  precisely  the  nonexistence  of  an  appropriate  front  that 
results  in  the  phenomenon  of  division. 

We  shall  be  examining  various  ways  to  relax  our  previous  assumptions,  and  shall  ash 
whether  the  above  necessary  condition  continues  to  hold.  Throughout,  we  retain  (1)  with 
a  =  e  as  the  fundamental  equation.  The  presence  of  the  factor  e  <<  1  with  ut  makes 
u  a  fast  variable.  The  role  of  this  parameter  in  the  foregoing  analysis  was  to  make  tne 
u-profile  respond  Quickly  to  changes  in  the  function  f.  Thus,  when  v  reaches  the 
critical  value  v_^,  this  quick  response  initiates  the  front  division  before  v  has  a 
chance  to  attain  the  other  critical  value  v,  after  which  time  a  front  reversal  w"uld 
occur,  hut  no  division.  Of  course  other  mechanisms  nay  delay  v’s  reaching  that  seven  i 
point.  Calculations  reported  in  [3]  indicate  division  may  indeed  occur  even  if  €  -  1. 
:i't  the  smallness  of  c  enhances  the  division  process,  and  we  shall  retain  that  assu-rtio 
•■hroughout  the  following. 
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The  parameter  e2  as  diffusion  coefficient  in  ( 1 )  can  always  be  obtained  by  proper  x- 
scaling,  and  we  retain  that  also.  The  only  real  effect  of  such  scaling  may  he  in  reaard  ", 
the  initial  data  and  Stage  1;  more  comments  on  this  will  be  given  in  the  next  section. 

(X)  Relaxation  of  the  symmetry  in  (3).  in  place  of  (3),  assume  (again  for 
simplicity) 


(  1  ♦ 
g(u)  -  ^  t 


s 


u  >  u0 
u  <  u0 


113) 


for  some  constant  s  >  0.  The  ideas  here  are  immediately  extendable  to  the  case  of  more 
general  discontinuous  function  g(u,v).  This  change  does  not  affect  Stage  1.  The  effect 
here  is  to  cause  the  front,  previously  fixed  during  Stage  2  at  x  «  0,  to  move.  Ir.  fact 
if  it  were  not  to  move  for  t  <  t0,  then  at  that  time  the  last  term  in  ( 1  )  could  be 
expressed  as 

f(  )  .  (f(u,v° '  V  '  u  4 

U  \  f(U'V0  *  t0(1+S)  ,  U  Ug  , 


,Vvo  +  V1+s) 

This  function  satisfies  Jh  (y  _  t  j  f(u)du  <  0,  hence  h_  is  the  dominant  state,  and 

the  front  will  be  moving  to  the  right.  So  it  cannot  be  that  the  front  was  stationary;  the 
effect  of  the  term  s  is  to  cause  it  to  travel  to  the  right. 

Let  x  ■  y(t)  >  0  be  the  position  of  the  front  (where  u  =  uQ).  A  further 
consequence  of  its  movement  is  that  v  is  no  longer  discontinuous.  In  fact  since  the 
motion  is  to  the  right,  the  value  of  v  at  the  front's  position  during  Stage  2  is 

v( y( t) , t)  =  Vq  +  t( 1+s)  , 

whereas  v(0,t)  *  Vq  -  t  for  the  same  reason  as  before.  Between  the  two  values  r>  a-  i 
y(t),  v  is  a  continuous  monotone  function  which  could  be  calculated.  This  continuous 
function  is  larger  than  the  discontinuous  function  equal  to  vn  -  t  for  x  <  y(t)  an' 

i 
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v0  ♦  t  (  1  ♦  s  ’  f  x  •  v‘  r  .  V  w  ^  if-  a  n1  {♦•iisir.ij  of  v,  ar  •*  *  •  *?*n* 

velocity  13  a  ifecreasi rm  ?  uo  :*■  ion  >f  f  ‘see  '*!,  Appendix',  so  t-hP  er  *•  -f  f  *Te  '»»'•*  l 
is  to  further  increase  the  \«e  i  or*  i  •  v  ,  -t-Mi  i  on  r 

This  acceleration  s  ^  w •  rat  t*e  asymnet  ry  ha  »  a  lest  ah  1 1  i  ?  .  ruj  ef^ec*:  >n  t  t.e  fror  .  ?*'*■ 

front  ’  s  velocity  dur*o.’  :*-a-n  ’  will  have  magnitude  h  f  s  j  '  s  *  *  ’  . 

We  must  asV  how  fai  ’(i*  ^ront  wi’’.  r--»v/*  dur !  n  »  Stage  2.  Since  this  stage  :s  or 
duration  ft),  t  •>*-.  1 1  at  .'••■'ce  w.ll  he  fs^.  As  «***en  from  the  above  argument,  this 
Interval  in  x  he*  •**«  firs?  ar.  1  final  positions  n*  the  front  in  S*-age  2  .s  wv'erc-  u 

of  the  spatial  van  iti:n  .r.  v  is  c '•ncer.tr  ate  1.  So  the  gradient  v  -  h(i),  t  sari 
the  condition  that  the  variation  he  0(1 l  in  a  front's  width,  however,  we  wust  have 

I v  |  >  0(1'.  This  : ,  n j  possible  if  s  *  Of  ), 
x 

We  con?lu?e  that  inder  Mi'.  front-  divisi.>"  .s  onlv  possible  if  s  <  h  <  \ .  i.  *■  ’•  »  s 

condition  is  fulfilled,  divisio-  .wur  the  firs*-  tine,  hut  ^f  course  subseouer* 

divisions  will  become  nor  -  uni  i «>*1  v. 

r  I X )  Relaxation  of  the  symmetry  ■  4  )  jvi_  f.  In  t^is  case,  even  thnuev  M>  ic 
retained,  it  will  no  Lon  ler  necessarily  be  true  that  the  front's  velocity  is  zen  durin 
Stage  2.  For  example  .consider  a  n.il  Inline  as  follows: 

-n 


If  the  front  were  artificially  held  fixed  for  0  <  t  <  tQ,  the  asymmetry  here  would 
try  to  push  the  front  forward,  because  /  f  du  over  the  discontinuous  line  shown  for 
t  -  t0  will  be  negative.  This  forces  the  front  to  move,  ha  before,  it  can  be  seen  that 
the  front's  movement  further  causes  acceleration,  so  the  asymmetry  is  a  destabilizing 
influence.  Let  s  be  some  measure  of  asymmetry  in  ti  for  example, 

.  u  h+(vQ  +  u)  - 

*  '  Hm,  (vQ  -  w)f(u'V0  "  “,<,U  *  /u0  f‘U,V0  +  U)dU|  ' 

the  maximum  being  over  all  u  for  which  h^tv^  *  u )  *re  both  defined. 

Then  the  argument  in  Case  (I)  can  be  repeated  to  show  that  front  division  is  only 
possible  if  a  <  0(c). 
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The  fact  that  v  is  not  constant  initially  as  in  (6)  will  have  little  effect  at  this 


stage.  What  is  more  important  is  the  place  (u.v)  where  r  crosses  the  nullcurve.  It 
was  assumed  previously  that  (u.v)  -  (Uq.Vq):  now  we  assume  this  is  not  so.  For  the  sake 
of  definiteness,  suppose  v  >  vQ  and  that  <|i  is  a  decreasing  function  of  x  'as 
indicated  in  Figure  6).  Since  h_  dominates  at  v  -  v,  we  have  c(v)  >  0,  so  at  the 
outset  in  Stage  2,  the  front  moves  to  the  right  at  a  nonzero  speed. 

Let  us  see  what  happens  in  the  simple  case  that  v(x)  is  linear,  so  ','(x)  »  -r  < 

Let  y(t),  as  before,  be  the  position  of  the  front.  Ahead  of  it,  v(x,t)  »  v'x)  ♦  t, 
x  >  y(t>,  hence  v+(t)  s  v(y(t),t)  satisfies 

**  •  ♦,(y)  y  ^  i  •  **  c  ♦  i  .  ns) 

Because  the  front  is  in  motion,  as  before,  the  function  v  will  not  be  discontinuous  at 
the  front,  and  the  velocity  c  will  be  the  trigger  velocity  associated  with  v  -  v*(t>. 
From  (15),  we  now  have 


v+  «  -me  ( v+ )  +  1  =  H  ( v+  )  .  '16) 

If  there  exists  a  value  v*  such  that  c(v*)  »  -,  it  will  be  a  rest  point  for  (16),  anl 

m 

since  c  is  an  Increasing  function  of  v,  H'(v+)  <  0,  implying  that  v*  is  a  global 
attractor.  Of  course  v*  >  Vg,  since  c(v*)  >  0.  Thus,  the  value  v  «  v0  is  never 
achieved  at  the  front,  and  the  front  does  not  divide  in  this  case. 

on  the  other  hand  if  -  is  so  large  that  there  is  no  trigger  velocity  ecrual  to  it, 

HI 

then  H(v)  >  0  for  all  v  e  <vQ,v).  A  little  further  analysis,  which  we  shall  not  1e*  i . 
then  shows  that  v+  grows  till  it  reaches  v,  at  which  point  the  front  accelerates  an 
becomes  a  "phase  front"  (5! ,  v+  remaining  at  v. 

Now  suppose  tt>*  is  not  constant,  but  rather  1'  ♦  constant  as  x  *  ®,  so 
Again,  it  is  clear  that  v+  -  y'c  +  1  >  0,  and  v+  increases  to  v,  wher»  it  «tay«. 
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All  this  was  under  *  he  asranjMor.  that  y '  <  On  the  other  hand  it  'Jd  >  0,  (15; 

shows  that  :.c  rest  point  exii/s,  an  i  v+  increast-r  for  a  fire.  It  never  reaches  vf 
however,  because  v  is  *nv.*ri»asin-i  a^  a  fur<*tjr  <jf  x,  and  before  v+  can  reach  v, 
therr  will  he  a  value  ..f  x  >  y(t)  a?  whi«*h  v  =  v.  At  such  a  point,  a  dovn^umr  wave 
front  (as  shown  in  a  similar  setting  in  [4, S’  )  is  .unerated.  It  mover:  to  the  left  and 
eventually  annihilates  t;ie  original  fr  >nt.  Again,  no  front  division  occurs. 

Up  to  now,  it  has  bee*  aar. omed  that,  the  characteristic  spatial  scale  of  the  initial 
data  is  :  uaf  associa-**  i  with  ♦  h»-  widtn  of  the  front;  so  <r  and  ^  are  well-behaved 
functions  of  -  -  ' .  This,  -or  example,  would  come  about  if,  in  the  original  equation  (1), 
a  *  i  and  »  ar  d  .  aie  functions  c*  k.  Now  suppose  the  initial  data  have  natural  scale 
?  when  t  *  '  : 


.  =  *•(  x/f  i ,  .  *  v<x/e>  . 

The  possible  effects  when  :  f  -  are  setrn  mainly  in  Staae  1.  If  p  <<  t,  there  is 

no  effect  at  all;  the  initial  function  will  vary  more  sharply  than  the  front  profile,  and 

equilibrat ior.  to  the  front's  profile  will  still  occur  in  a  time  interval  O'nl.  (This  was 
true,  in  fact,  in  M21,  wher*  tre  initial  data  was  discontinuous,  corresponding  to 

2  +  0.)  on  twe  oth*r  Mr.  i,  -he  initial  sharp  variation  in  v  may  be  such  as  to  maVo  :i 

front  division  imminent,  .er-en  f  course,  or.  the  initial  data.  T^e  irguments  above 

sv'ow  t^af  f  >r  -hi«-  *'  happ*  r ,  -•••  ••  i  r  ial  v  distrihitiog  r>oul  i  he  such  a~  to  -naVe 
vha  :  a~t«r  ;  st  i  *  fron*  Vf!"'-;*  ?'*•'■»  a-r  r  ~  <  i^ar^l  y  sm  T'- .  s  ,  again  ,  '  n  a  special 

con  ii  -  ion . 

If  d  >  >  ,  howftv  v  ,  M'  -*<  pt-i  i^av*'-  of  -  i  -  j »-  1 1  data  ~r  ,ss.»s- he 

j  s -*ur»  !  i ’»7  ‘’•Mr.'-1,  in  *  »  *"  '  s.  ,  i»,  -  ho  «*on  :  l  ll  rat  i  on  in  c  -  a  a ''  1  will  bo  dolaved* 

This  is  Vf*  •'  t  : :  no  o-  >  -n  -ra-:a'  .  ’  .» n  lo*'  :  if'*!  with  the  *rnr;t*  ••tmi  1  i *  r  :  urn  prorijo, 

u  Will  be  near  )  y  cor  v  in-  i*  1  *  *.»'  -  .  '  i r  ■<  ur)  \  >  a  r*.'^  ■’fate  My-  u  in  f  r<  ) 
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(though  an  unstable  one),  it  will  take  a  long  time  for  u  to  leave  a  small  neighborhood 
of  uq ,  resulting  in  the  mentioned  delay.  This  delay  will  not  enhance  the  process  of  the 
front's  division,  and  may  in  fact  prevent  it. 

Finally,  some  comments  should  be  made  about  the  special  form  of  (2),  namely  that  v 
does  not  diffuse  and  that  g  is  a  discontinuous  function  of  u.  Both  of  these  provisions 
were  for  the  purpose  of  making  the  analysis  simpler  and  to  provide  for  v' s  sharp 

variation.  The  sharp  variation  will  still  be  the  case  if  v  is  allowed  to  diffuse  with  a 

2 

small  diffusion  coefficient,  <  0(£  ),  or  if  g  is  allowed  to  be  continuous,  but  with  a 
large  enough  gradient  near  u  «  ufl.  Me  shall  not  pursue  the  details.  In  [2]  and  [3],  in 
fact,  calculations  were  performed  with  a  continuous  g.  If  v  diffuses  by  a  relatively 
large  amount,  then  it  will  not  be  able  to  achieve  the  variation  0(1)  within  the  confines 
of  a  narrow  front,  and  so  will  not  support  a  division. 
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r> .  Discussion 

The  phenomena  associated  with  Yakhno's  mechanism  are  striking  and  suggestive,  but 
rather  special,  as  they  depend  on  a  great  deal  of  symmetry  in  the  basic  equations.  A 
disruption  of  this  symmetry  will  generally  destroy  even  the  qualitative  features  of  the 
phenomenon.  In  order  to  serve  as  a  realistic  model  of  leading  centers  in  biological 
contexts,  the  basic  model  will  have  to  be  supplemented  by  other  mechanisms,  possibly 


through  additional  reaction-diffusion  equations,  which  will  serve  to  restore  structural 
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